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PLASTIC  DEFORMATION  OF  A  THIN  CIRCULAR  PLATE  UNDER  PRESSURE 


ABSTRACT 

In  this  report  a  numerical  solution  is  given  of  a  set  of  equations 
consisting  essentially  of  three  plasticity  laws,  two  strain-displacement  laws, 
and  two  equilibrium  laws,  which  describe  the  action  of  a  clamped,  thin,  cir¬ 
cular  plate  as  it  yields  plastically  when  pressure  is  applied  to  one  side. 

The  stresses,  strains,  thickness  variation,  and  deflections  for  any  thin  cir¬ 
cular  plate  of  a  given  material  may  be  computed  by  the  numerical  integration 
of  the  equilibrium  conditions,  the  geometric  conditions  relating  displace¬ 
ments  and  strains,  and  the  stress-strain  laws.  The  solution  may  be  reduced 
to  the  solution  of  a  second-order  differential  equation  with  the  radial  dis¬ 
tance  r  as  independent  variable.  The  solution  depends  upon  an  experimentally 
determined  function,  r(y),  which  describes  the  stress-strain  properties  of 
the  material,  and  upon  three  parameters,  the  pressure  p,  the  original  thick¬ 
ness  h0,  and  the  radius  a  of  the  clamping  ring.  It  is  found  that  for  a  given 
material,  a  family  of  curves  with  pa/h0  as  a  parameter  serves  to  predict  the 
solution  for  any  thin  circular  plate  of  the  same  material. 

This  analysis  has  been  carried  out  for  a  particular  function  r(y) 
based  on  results  of  a  tensile  test  made  on  a  specimen  of  medium  steel.  Graphs 
of  theoretically  and  experimentally  determined  values  of  deflection,  radial 
and  circumferential  strains,  radial  and  circumferential  stresses,  and  thick¬ 
ness  corresponding  to  various  pressures  are  presented  which  apply  to  all 
plates  made  of  the  same  steel  as  this  specimen. 

INTRODUCTION 

The  plastic  action  of  a  circular  plate  under  static  pressures  great 
enough  to  cause  rupture  is  of  Interest  to  the  David  Taylor  Model  Basin  in  its 
program  of  study  of  damage  to  ships  caused  by  underwater  explosions.  The 
clamped  circular  plate  is  a  simple  model  regarded  as  analogous  in  some  re¬ 
spects  to  one  of  the  plate  panels  in  a  ship.  The  static  pressure  applied  to 
the  circular  plate  is  in  the  ship  replaced  by  the  pressure  due  to  an  under¬ 
water  explosion. 

When  sufficient  pressure  is  applied  to  one  side  of  a  circular  steel 
plate  clamped  around  its  edge,  the  plate  is  deformed  Into  an  approximately 
spherical  dome  or  cap.  If  the  clamping  ring  has  a  curved  edge  as  in  Figure 
1 ,  the  rim  of  the  cap  is  tangent  to  the  curved  edge .  As  the  plate  deforms 
the  metal  flows  so  that  the  plate  thins  markedly  at  its  center  and  only 
slightly  near  the  rim.  A  point  P  on  the  initially  flat  plate  moves  along 
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Figure  .1  -  Diagram  Showing  a  Thin  Circular  Plate 
as  It  Deforms  Plastically  under  Pressure 

» 

a  curve  such  as  PP'  to  P'.  For  a  given  deflection  of  the  diaphragm  the  ra¬ 
dial  displacement  u  is  zero  at  the  center  by  reason  of  symmetry,  reaches  a 
maximum  somewhere  between  the  center  and  edge  of  the  plate,  and  drops  to 
zero  again  at  the  rim. 

The  behavior  of  a  circular  plate  in  the  elastic  range  has  been 
considered  by  Hencky  (1)*  and  by  S.  Way  (2).  In  Hencky's  solution  bending 
stresses  are  neglected.  Way  has  taken  both  bending  and  stretching  into  ac¬ 
count.  Since  the  present  analysis  is  concerned  with  thin  plates,  the  bending 
stresses  are  neglected  here  also.  This  analysis  differs  from  Hencky's  in 
that  it  is  concerned  primarily  with  the  plastic  range  up  to  rupture.  The  re¬ 
sults  presented  here  are  valid  for  values  of  z0/a  ranging  from  about  0.03  up 
to  about  0.30,  where  z0  is  the  central  deflection.  For  values  of  z0/a  cor¬ 
responding  to  elastic  deflection,  that  is,  when  z0/a  is  less  than  0.03,  the 
results  nevertheless  have  some  validity  since,  as  will  be  shown,  the  plastic¬ 
ity  laws  chosen  are  identical,  for  small  strains,  with  the  elasticity  laws 
provided  Poisson's  ratio  Is  taken  as  0.5. 

The  purpose  of  the  present  report  is  to  show  how  certain  plastic¬ 
ity  laws,  three  in  number,  which  involve  an  experimentally-obtained  octahedral 
stress-strain  curve  r(y),  may  be  used  to  compute  stresses,  strains,  and  dis¬ 
placements  in  a  thin  circular  plate  when  merely  the  dimensions  of  the  plate, 
the  pressure  applied,  and  the  true  tensile  stress-strain  curve  up  to  rupture 
of  a  specimen  of  the  material  are  given.  The  present  problem  has  its  paral¬ 
lel  In  elasticity  theory,  but  the  methods  employed  in  this  report  are  almost 
necessarily  graphical  and  numerical  due  to  the  nonlinear  character  of  the 
stress-strain  laws  of  plasticity.  Emphasis  will  be  placed  on  finding  a  fea¬ 
sible  approximate  method  to  solve  the  general  plasticity  problem  of  the 
clamped,  thin,  circular  plate  rather  than  on  attaining  a  close  check  between 
experiment  and  theory.  Improvements  in  the  theory  which  would  yield  numer¬ 
ical  results  in  closer  agreement  with  experimental  results  are  suggested  at 
the  end  of  the  report. 


*  Numbers  in  parentheses  indicate  references  on  page'  30  of  this  report. 


BASIC  ASSUMPTIONS 

The  calculation  of  stresses,  strains,  and  displacements  involves 
solving  simultaneously  two  equations  which  express  the  conditions  of  equilib¬ 
rium,  two  geometrical  relationships  of  strain  and  displacement,  and  three 
plasticity  laws  which  relate  stress  with  strain  in  the  material.  In  the  theo¬ 
ry  the  plate  is  assumed  to  be  clamped  by  a  rigid  ring  with  a  sharp  inner  edge; 
bending  moments  are  assumed  negligible  since  the  plate  is  thin.  Consequently, 
the  surface  of  the  loaded  plate  meets  the  plane  of  the  rim  at  a  finite  angle. 
The  methods  and  formulas  used  in  the  calculations,  excepting  perhaps  the  plas¬ 
ticity  laws,  are  exact  but  for  the  omission  of  second-order  terms.  Hence  the 
validity  of  the  results  obtained  hinges  largely  on  the  validity  of  ®the  three 
plasticity  laws  which  are  assumed.  Of  these  laws,  one,  which  expresses  con¬ 
stancy  of  volume  of  the  material,  has  been  verified  experimentally  with  a 
high  degree  of  accuracy  for  certain  metals.  Agreement  with  experiment  is 
reasonably  good  for  the  other  two  laws. 

DEFINITION  OF  STRAIN 

Although  there  is  no  question  in  the  literature  as  to  the  defini¬ 
tion  of  stress  and  the  stress  tensor,  there  is  some  choice  in  the  definition 
of  strain.  The  conventional  strain  e  =  — -I— Lp ,  where  l  and  £0  are  the  final 
and  the  initial  lengths  respectively,  may  be  considered  basic  and  other  defi¬ 
nitions  of  strain  can  be  made  in  terms  of  it.  The  logarithmic  strain,  for 
example,  designated  as  e' ,  has  been  defined  by  Ludwik  (3)  as 

£ 

e'  ~  Jr+7  "  |0B(1  +  f) 

o 

where  e  is  the  conventional  strain.  Another  definition  based  on  the  quad¬ 
ratic  form  for  the  square  of  the  differential  of  arc -length,  one  which  guar¬ 
antees  the  existence  of  a  strain  tensor,  has  been  considered  by  Synge  and 
Chien  (4).  For  the  sake  of  simplicity,  and  because  the  equations  to  be  em¬ 
ployed  are  necessarily  approximate  themselves,  strain  will  be  taken  to  mean 
conventional  strain  in  this  report. 

STRESS-STRAIN  EQUATIONS  OF  PLASTICITY 

The  plasticity  laws  stated  in  the  following  were  given  by  Nadai  (5) 
(6);  they  are  based  on  the  work  of  a  number  of  writers,  among  whom  may  be 
mentioned  von  Mises  (7)  and  Mohr  (8).  White  (9)  has  applied  them  to  a  theo¬ 
retical  investigation  of  circular  diaphragms  under  pressure.  They  have  been 
tested  experimentally  by  Schmidt  (10),  Taylor  and  Qulnney,  Ros  and  Eichinger, 
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Nadai  and  Lode,  and  others.  Nadal  (6)  {11 )  has  written  accounts  of,  and 
given  references  for  the  theoretical  and  experimental  work  of  these  and  other 
writers.  Davis  (12)  (13)  has  made  a  series  of  tests  on  cylindrical  tubes  to 
verify  these  plasticity  laws  for  copper  and  medium  steel. 

The  stress-strain  relations  of  plasticity  which  are  assumed  in  this 

report  are 


+  e2  +  ««  "  0  [la] 

£u=_£i  »  =  £3_z_ii  =  1  [1b] 

<}\  <t2  <?2  &3  <7]:  X 

*  r  =  r(y)  [  1 C  ] 

where 

r  =  ~S  i^1  ”  +  ~  az)2  +  ^  ~  ffJ2]2  [2a] 

y  =  ~  *z)2  +  («2  ~  *3)2  +  (^3  _  fi)2  2  [2b] 

and,  as  will  be  shown. 


A  -  2^  [3] 

In  these  equations  elt  e2,  and  e3  are  the  principal  strains  and  ax,  <r2,  and 
aa  are  the  principal  stresses.  The  quantities  y  and  r  have  been  shown  by 
Nadai  (11)  to  be,  respectively,  the  shear  strain  and  the  shear  stress  in  an 
octahedral  plane,  a  plane  whose  normal  makes  equal  angles  with  the  three  prin¬ 
cipal  directions  of  strain  or  stress. 

Equation  [la]  implies  that,  to  a  first  order,  volumes  remain  con¬ 
stant.  Equation  [lb]  is  usually  associated  with  the  assumption  that  at  any 
one  point  the  principal  axes  of  stress  coincide  with  the  principal  axes  of 
strain.  It  may  be  shown  that  the  latter  condition  is  satisfied  for  a  thin 
circular  plate  under  pressure.  The  quantity  X  depends  on  the  material  as 
well  as  on  the  state  of  stress  and  strain  at  each  point.  Equation  [1c]  im¬ 
plies  that  v( y)  is  the  same  function  of  y  for  all  states  of  stress  and  strain 
In  a  material.  Thus,  if  y  has  the  same  numerical  value  in  two  states,  then 
by  this  formula  r  will  have  the  same  value  in  the  two  states.  Hence  r(y)  may 
be  considered  to  characterize  the  material  with  respect  to  its  stress-strain 
properties.  It  can  be  determined  by  a  tensile  test  on  the  material,  or  by 
tests  on  plates  deformed  under  pressure,  using  formulas  given  later  in  this 
report . 

It  may  be  noted  that,  since  e3  *  -  («j  +  «2), 
y  =  2 J/|-  \/*\  +  ej  f2  +  «22 


[4] 
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Equations  [1]  are  three  Independent  relations  which  connect  the 
six  variables  clt  e2,  e3  and  alt  a2l  <r3.  They  correspond  to  the  three  stress- 
strain  laws  of  elasticity,  and,  in  fact,  may  be  put  into  a  form  which  corre¬ 
sponds  closely  to  these  laws.  Combining  Equation  [lb]  with  Equations  [2a]  and 
[2b],  it  is  easily  seen  that 


k  = 


and  as  a  consequence  of  Equation  [1c]  the  ratio  k  is  seen  to  be  a  definite 
function  of  y  also.  The  principal  shear  stresses  may  therefore  be  expressed 
in  terms  of  shear  strains  by  the  formulas 


a\  ~  <72 
2 

a2  ~  g3 

2 


r(y) 

y 


(«1 


\ 


*2) 


g3  a\ 


2 


r(y) 

y 


(f3 


«l) 


[5] 


These  equations  are  equivalent  to 


y 

= 

y 

=  fM. 
y 


+  x 

e2  + 

C3  +  L 


[6] 


where  D  =  ■^(a1  +  a2  +  <r3 )  is  the  mean  normal  stress  as  may  be  seen  by  equat¬ 
ing  the  sum  of  the  three  left-hand  members  of  Equations  [6]  to  the  sum  of  the 
three  right-hand  members  of  these  equations  and  making  use  of  Equation  [la]. 
The  physical  Implication  is  that  a  hydrostatic  stress  may  be  added  to  exist¬ 
ing  stresses  without  changing  the  strains. 

Equations  [6]  solved  for  strains  in  terms  of  stresses  give 


*3) 


> 


11) 
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where  y(r)  is  the  function  which  is  inverse  to  r(y) .  These  expressions  are 
of  the  same  form  as  the  laws  of  elasticity;  the  factor  y{r)/J>r  replaces  1  /E 
and  the  factor  1/2  replaces  v,  where  v  denotes  Poisson's  ratio. 

DETERMINATION  OP  r(y)  PROM  A  TENSILE  TEST 

The  function  r{y)  may  be  expressed  in  terms  of  quantities  measured 
in  a  tensile  test  of  the  material.  This  is  true  because  in  such  a  test,  if 
the  three  principal  stresses  are  denoted  by 

=  °L  a2  =  0.  *8  =  0 


Figure  2  -  A  Load-Strain  Curve  for  a  Medium-Steel  Tensile  Specimen 

This  load-strain  curve  was  obtained  from  a  round  medium-steel  specimen,  0.505  inch  in 
diameter,  made  from  the  same  sheet  of  medium  steel  as  the  experimental  plate  described 
in  this  report.  This  curve  is  considered  to  "describe11  the  material;  it  was  used  as 
the  basis  for  all  computations  for  the  theoretical  prediction  of  the  action  of  a  thin 
plate  of  this  material. 
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and  the  three  principal  strains  by 


—  clj 

1  1 
£t  ~  ~  ~2  6i>  e3  ~  ~  2  61 

then,  from  Equations  [2a]  and 

[2b] 

£|« 

II 

h 

[8a] 

vy 

*2 

ii 

[8b] 

Since  the  load  L  and  the  strain  ex  are  the  quantities  usually 
measured  in  a  tensile  test,  it  is  useful  to  write  the  equation  for  ax  in 


Figure  3  -  Octahedral  Shear  Stress  Plotted  against  Octahedral  Shear  Strain 

The  calculations  of  this  report  were  based  on  the  r(y)  curve  shown  here.  The  load-strain  curve 
from  which  t( y)  was  computed  is  shown  in  Figure  2.  The  r(y)  curve  calculated  from  measurements 
at  the  center  of  an  experimental  plate  of  the  same  material,  medium  steel,  is  shown  for  compari¬ 
son.  The  difference  between  the  two  curves  accounts  in  large  part  for  the  difference  between 
computed  and  experimental  values  of  deflection  in  Figure  13  on  page  28. 

The  quantities  r  and  y  for  the  tensile  test  were  computed  by  Equations  [8]  and  [9],  For  the 
plate  test  they  were  computed  by  Equations  [11  ],  [16],  [17],  [18],  and  [19]  by  use  of  measured 
values  of  pressure,  deflection,  and  strains.  For  strains,, in  the  region  of  necking  of  a  tensile 
specimen,  r(y)  should  be  determined  from  measurements  of  instantaneous  load  and  minimum  cross- 
sectional  area  A  of  the  specimen.  In  this  region  should  be  taken  as  L/A  and  as  (A0  -  A) /A. 
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terras  of  them.  In  a  tensile  test  the  cross-sectional  area  is 

A  =  i40(l  +  e2 )  (1  +  ^3) 

where  A0  is  the  original  cross-sectional  area.  It  follows  that 


0,  =  —  = 


L_ 

A 


L_ 

A0 


1 


M'«.) 


Therefore,  approximately, 


<ti  = 


f  (1  +  Cl) 

Ao 


[9] 


The  curve  r(y)  of  Figure  3>  which  shows  the  relation  between  octa¬ 
hedral  shear  stress  and  octahedral  shear  strain,  was  computed  from  the  load- 
strain  curve  L{ex)  of  Figure  2  by  the  use  of  Equations  [8]  and  [9].  For  com¬ 
parison,  the  r(y)  curve  determined  by  measurements  on  an  experimental  plate 
of  the  same  material  under  pressure  is  included  in  Figure  3.  The  extent  to 
which  the  two  curves  agree  is  a  measure  of  the  validity  of  the  plasticity 
laws  assumed  in  this  report.  The  difference  between  the  two  curves  accounts 
for  a  large  portion  of  the  difference  between  the  experimental  and  the  theo¬ 
retical  results  presented  herein. 


Figure  4  -  The  Principal  Stresses 
and  Strains  in  a  Circular 
Plate  under  Pressure 


STRESS-STRAIN  LAWS  FOR  A  THIN  PLATE 

As  a  thin  circular  plate  clamped 
at  the  rim  deflects  under  pressure,  sym¬ 
metry  requires  that  one  of  the  principal 
v stresses  be  the  stress  ar  In  the  radial 
direction  and  another  be  the  stress  ae in 
the  circumferential  direction,  as  in  Fig¬ 
ure  4.  Similarly  one  of  the  principal 
strains  is  er,  the  strain  in  the  radial 
direction,  and  another  is  e9,  the  strain 
in  the  circumferential  direction.  If  ox 
and  <s2  of  Equations  [6]  are  interpreted 
to  be  ar  and  a6  respectively,  then  az  is 
the  stress  in  the  direction  perpendicular 
to  the  plate  at  the  point.  This  stress  is 
of  the  order  of  the  pressure  on  the  plate 
and  hence  is  of  relatively  negligible 
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,  magnitude.  If  a3  is  set  equal  to  zero  and  the  relation  «3  ■  -(ej  +  e2)  is 
used  in  Equation  [5],  the  result  will  be 

T  (v) 

<rr  =  2-yZL(2€r  +  «9)  [10a] 

(2««  +  «r)  [10b] 

where  now,  by  Equation  [4], 

y  -  2(/|^V€7~+17e7^r^?  [11] 

Also,  for  biaxial  stresses  Equation  [2a]  implies  that 

r  =  -y- 

The  foregoing  relations  may  be  called  the  plastic  stress-strain  equations 
for  a  plate  under'  biaxial  stress. 

STRAINS  IN  TERMS  OP  DISPLACEMENTS  AND  DEFLECTIONS 

Suppose  a  flat  circular  plate  is  deflected  to  a  radially  symmetrical 
shape  as  shown  in  Figure  1 .  The  axis  of  symmetry  may  be  taken  as  the  z-axis 


Figure  1  -  Diagram  Showing  a  Thin  Circular  Plate 
as  It  Deforms  Plastically  under  Pressure 

in  a  set  of  cylindrical  coordinates  r,  6,  z.  As  the  plate  deforms,  a  point  P 
on  the  plate  moves  along  a  path  such  as  PP\  Let  u  denote  the  radial  dis¬ 
placement  for  the  point  P\  If  *  is  the  deflection  of  this  point  then  the 
point  P  of  coordinates  r0,  0O,  0  on  the  undeformed  plate  has  the  final  coor¬ 
dinates  r,  0,  z.  By  symmetry  e  »  0O  and  the  square  of  the  arc  length 

ds0z  -  dr*  4-  r0sd$ 02 

of  an  element  on  the  undeformed  plate  has  as  its  final  arc  length 


ds2  —  dr 2  +  r2d$2  +  di2 
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The  strain  for  this  element  is  therefore 


ds  —  ds0  _  ds_  _  j  _  /dr2  +  r2d$L  +  dz‘ i 
ds0  ~  ds0  ~  y  dr0 2  +  rQ2d602 


[12] 


If  d$  =  d90  is  set  equal  to  zero  in  this  expression,  the  radial  strain  is 
found  to  be 


or, 


dr1  -f  dz‘ 
drn 2 


-  1 


Therefore,  since  r0  =  r  -  u, 


[13] 


If  du/dr  and  dz/dr  are  small  quantities.  Equation  [13]  may  be  Witten  in  the 
approximate  form 


_  du_  1  (dz\2 

€r  dr  2 \dr> 


[14] 


By  setting  dr0  =0,  and  dr  *=  dr0  +  du  -  du  in  Equation  [12]  the 
circumferential  strain  is  found  to  be 


du2  +  r2d0 2  +  dz2 


r02  de02 


-  1 


Therefore,  since  6  «  90i 


Because  of  radial  symmetry,  du/d$  -  0  and  dz/d8  *  0,  Hence 


U  -  “  -  1 


n 


<r» 


t 


or 


’■o 


In  terms  of  the  final  radius  r, 

•  u 


Since  u  is  small  relative  to  r,  we  write  for  the  purposes  of  this  report 
simply 


= 


u_ 

r 


[15] 


Equations  [ 1 4 ]  and  [15]  for  strains  in  terms  of  displacements  are 
the  familiar  ones  for  the  circular  plate  in  elasticity  theory. 


THICKNESS  IN  TERMS  OP  STRAINS 

Let  et  denote  the  strain  in  the  direction  normal  to  the  surface  of 
the  plate.  Then  by  Equation  [la] 


€z  —  —  («r  +  *i) 

If  h0  is  the  initial  thickness  and  h  the  thickness  of  the  plate  when  de¬ 
formed,  then,  by  definition. 


It  follows  that 


h  =  h0  (1  —  'er  -  et) 


[16] 


EQUILIBRIUM  CONDITIONS 

Suppose  that  under  a  pressure  p  the  deflection  function  for  a  cir¬ 
cular  plate  is  z(r).  The  force  resulting  from  the  pressure  on  a  cap  of  radius 
r,  Figure  5,  in  the  direction  of  the  z- 
axis  is  7 rr2p.  This  force  is  balanced 
by  a  vertical  force*  due  to  the  radial 
tension  tr  in  the  plate: 


2  nr 


dz 

dr 


Figure  5  -  Forces  on  a  Spherical  Cap 


*  The  sign  of  the  square  root  which  appears  in  this  expression  is  chosen  so  that  the  expression  is 
positive. 


12 


When  these  forces  on  the  cap  are  equated, 
there  results 


n  r2p 


—  2nr 


dz 

dr 


vM 


dz_\2 
dr  / 


or 


tr  =  —P 


2 


[17] 


dr 


Figure  6  -  Forces  on  an  Element 
of  a  Spherical  Cap 


Consider  now  the  element  defined  by  A9  and 
As  as  in  Figure  6.  The  resultant  force  on 
the  element  due  to  the  circumferential  ten¬ 
sions  te  at  each  end,  as  shown  in  Figure  6, 
is  approximately  teA9As.  This  force  acts 
radially  Inward  in  a  plane  parallel  to  the 
original  plane  of  the  plate.  The  component  of  this  force  in  the  direction  of 
the  radial  tension  is  approximately  tgA0  As  jf  ~  tg  A9  Ar.  The  total  force  on 
the  element,  due  to  the  radial  tension  tr>  is  approximately  (r  +  Ar)  AO  ( tr  + 
Atr)  -  rA0tr.  No  component  due  to  pressure  exists  in  the  tr  direction.  Hence, 
if  the  forces  in  the  tr  direction  are  equated,  there  results  the  approximate 
equation 

tgA9Ar  =  (r  4-  Ar)A9(tr  +  Atr)  ~  rA9tr 


or 


to  — 


(r  +  Ar){tr  +  Atr)  -  rtr 


°9  ~  Ar 

If  we  let  Ar  approach  zero,  we  obtain  the  exact  equation 

U  =  Jf{rtr)  "  tr  +  r  TZ^t) 


dr 


[18] 


or 


~  -7TP 


dzz 
dr * 


©> Ml 


Thus,  both  tr  and  t9  may  be  expressed  in  terms  of  the  deflection 
function  z[r)  and  the  pressure  p.  The  stresses  og  and  af  may  be  expressed  in 


l 


13 


terms  of  tr,  te,  and  the  thickness  h  as  follows, 


tr 

*'  =  T 

[19a] 

it 

£ 

[19b] 

By  combining  Equations  [ 1 9 ]  with  Equations  [14],  [15],  [16], 
it  is  seen  that,  if  the  displacements  u  and  z  of  every  point 
the  pressure  p  are  known,  the  stresses  may  be  computed. 

[17],  and  [18], 
on  the  plate  and 

SOLUTION  OP  EQUATIONS  FOR  PLASTIC  DEFORMATION  OP  A  CIRCULAR  PLATE 

Prom  the  equations  derived  in  the  foregoing,  the  following  Inde¬ 
pendent  set  of  equations  may  be  chosen  to  describe  the  action  of  the  thin 
circular  plate  unaer  pressure: 

*>  - -Hr +  (£)‘ 

dr 

[17] 

te  -  ir  +  r  £p(tr) 

[18] 

<rr  =  2  {2er  +  €#) 

[10a] 

Og  2  y  (2  € g  +  €  r) 

[10b] 

du  .  1  (dz\z 

6r  “  dr  2  l ir) 

[14] 

u 

»  r 

[15] 

h  **  Kq  (l  •—  €f  £$) 

[16] 

tr 

r  T 

[19a] 

$ 

11 

[19b] 

>V  +  *«,+  «•) 

[11] 

14 


These  equations  are  ten  Independent  relations  in  the  ten  variables 
t  r ,  t$,  z,  u ,  <xr,  <?e ,  «r,  €g,  y,  and  h,  containing  the  arbitrary  parameters  p 
and  h0  and  the  function  r(y) ,  which  describes  the  material  with  respect  to  its 
plastic  stress-strain  properties.  The  equations  will  now  be  transformed  into 
others  suited  to  solution  by  numerical  integration. 

Equations  [ 1 4 ]  and  [153  imply 


or, 


[20] 


By  solving  for  in  Equation  [17  ] ,  we  find  that 

dz  _  vr 

dr  ~  yuS  -  p2r2 

so  that  Equation  [20]  becomes 


l  de9 
+  r  -~ 
dr 


1  p2r2 

2  4  tr  —$rl 


[21] 


[22] 


Equations  [10]  and  [19]  imply 


tr  _  2tr  +  6y 

t9  2c9  er 


or,  by  Equation  [ 1 8 ] 


Thus 


tr( 2(g  4  €r)  — 


tr  +  r 


d(tr) 
dr  J 


( 2  f  r  +  e.) 


[23] 


(r<€,  -  ,r>  =  r  +  ',) 

or 

d(tr)  _  ~  Gr  tr  f  .  . 

dr  2cr  +  ee  r  [24J 

Prom  Equation  [10a]  and  Equations  [1 6 ]  and  [19]  there  follows 

tr  -  K  2  ^(2er  +  *«)(1  -  €r  -  e9)  [253 

Equations  [22],  [24],  [25],  and  [11],  are  four  equations  in  the  four  vari¬ 
ables  €r,  e9,  y ,  and  tr.  Inspection  shows  that  they  are  equivalent  to  a 


15 


*  * 


single  second-order  differential  equation  in  one  dependent  variable,  say  e9. 
There  will  consequently  be  two  arbitrary  constants  of  Integration.  These 
constants  vfill  be  determined  by  two  boundary  conditions,  for  example,  e9(a)  = 
0,  ee(-a)  =  0,  where  a  is  the  prescribed  radius  of  the  plate.  The  solutions 
of  these  equations  and  Equation  [21]  will  depend  not  only  on  the  radius  a,  but 
also  on  the  two  arbitrary  parameters  p  and  h0.  Fortunately,  the  triple  in¬ 
finity  of  solutions  resulting  may  easily  be  obtained  from  a  single  infinity  of 
solutions  depending  on  only  one  parameter  pa/h0  .  To  verify  this  statement 
it  is  necessary  merely  to  rewrite  Equations  [21],  [22],  [24],  [25],  and  [11] 
in  the  form 


1  tr  £g 

r_  h0  2er  +  e9 
a 


£  =  2  ^-(2er  4-  e,)(l  -  er  -  eg) 

n0  y 


y  -  2j/y"  /(er2  +  er€g  +  eg2) 

In  the  foregoing  equations  r/a  is  the  independent  variable,  and  the  dependent 
variables  are  now  z/a,  cr,  efl,  y,  and  tr/h0,  the  parameter  being  pa/h0  .  The 
two  boundary  conditions  are  that  e9  =  0  when  |  -  ±1.  Although  the  existence 
of  these  "proportionate  coordinates"  has  now  been  established,  it  is  con¬ 
venient  to  continue  the  analysis  with  the  original  coordinates.  Only  the 
graphs  of  the  solutions  need  be  given  in  proportionate  units. 

After  the  functions  er(r),  eB[r),  and  tr(r)  have  been  determined  by 
solving  Equations  [22],  [24],  [25],  and  [11],  the  deflection  function  z(r) 
will  be  found  by  integrating  Equation  [21 ]j  the  constant  of  integration  is 
fixed  by  the  boundary  condition  2=0  when  r  =  a.  The  quantities  such  as  t6 , 
<rt,  <tBs  w,  and  h  may  also  be  computed  by  Equations  [10],  [11],  [15],  [16], 
[17],  [18],  and  [19]. 
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It  is  convenient  to  plot  tr/h0  as  a  function  of  cr  and  «9,  using 
Equations  [11]  and  [25].  Figure  7  is  such  a  graph,  in  which  r(y)  is  taken  as 
in  Figure  5,  which  in  turn  was  based  on  the  tensile  test  curve  of  Figure  2. 
Figure  7  may  be  considered  to  describe  the  biaxial  stress-strain  properties 
of  a  thin  plate  of  the  material. 


NUMERICAL  SOLUTION  OF  EQUATIONS  [11],  [21],  [22],  [24],  [25] 
Equations  [24]  and  [22]  are  equivalent  to 


tr 


*0 


/■ 


tr  er _ €  % 

r  2er  +  e9 


dr 


[26a] 


.  [( 1  P2r 

e°  J  l  2  4 trz  —  p2r2 


[26b] 


where  tn  is  the  value  of  the  tension  tr  at  r  =  0  and  e0  is  the  value  of  the 
strain  e9  at  r  ■  0. 

Equations  [26]  are  solvable  numerically  by  iteration  of  the  two 
integrals,  by  the  use  of  Equations  [11]  and  [25],  as  follows:  Suppose  ap¬ 
proximations  for  tr  and  e9  in  terms  of  r  are  known.  Then  Equations  [11]  and 
[25],  or  the  graph  of  Figure  7,  give  er  in  terms  of  r.  Substituting  these 
approximations  of  tr,  e9,  and  cr  in  the  right-hand  side  of  Equations  [26]  and 
integrating,  we  obtain  new  approximations  for  tr  and  ee.  The  process  may 
then  be  repeated  indefinitely.  A  measure  of  the  accuracy  of  a  solution  is 
its  proximity  to  the  previous  one. 

When  sufficiently  accurate  evaluations  of  er  and  eg  have  been  se¬ 
cured,  Equations  [10]  may  be  used  to  find  ar  and  ae  . 

Finally,  the  deflection  function  z(r)  may  be  evaluated  by  an  inte¬ 
gration,  since  by  solving  for  z  in  Equation  [21],  there  results 


Pr 


where  z0  is  evaluated  so  that  z(<t)  *  0.  Then  z6  is  the  calculated  central 
deflection  of  the  plate. 


CHOICE  OF  PARAMETERS 

A  number  of  details  concerned  with  the  process  of  the  numerical 
solution  are  of  interest.  It  has  been  seen  that  there  are  but  three  independ¬ 
ent  parameters,  Hence  only  three  of  the  six  quantities  «0,  ta,  c0,  h0,  p,  and 
a  may  be  chosen  in  advance.  Thus,  for  example,  the  thickness  and  radius  of  a 
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Figure  7  ■  Tension  Due  to  Biaxial  Strain 

Let  and  e2  be  the  principal  strains  in  the  plane  of  a  thin  plate  and 
let  tj.  and  t2  be  the  corresponding  principal  tensions.  The  curves  shown 
here  are  curves  of  constant  fiA0  plotted  on  ex  and  t2. 

The  quantity  t iAa  corresponding  to  <x  and  t2  was  computed  by  the  formulas 

T ~  “  2  — ^-(2«1  +  f2)(l  —  £]  —  £2) 

Ho  y 

y  -  2 +  <i«2  +  <22 

The  function  r(y)  was  taken  to  be  the  curve  of  Figure  3  obtained  from  the 
tensile-test  curve  on  a  medium-steel  specimen  shown  in  Figure  2;  see  pages 
6  and  7  respectively. 

If  any  two  of  the  three  quantities  t\/h0,  elf  and  £S  are  given,  the  third 
quantity  may  be  found  by  interpolation  on  the  graph. 

This  graph,  where  we  have  taken  tr  =  tlr  t  f  =  tlt and  <$  =  *2,  was  used  to 
find  cr  in  the  theoretical  computations  from  tr  and  <$. 
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plate  may  be  given  and  the  deflection  and  strains  may  be  sought  as  functions 
of  the  pressure  applied.  Or,  in  other  instances,  the  radius  of  a  plate  and 
its  central  deflection  at  a  given  pressure  may  be  fixed  and  the  thickness  re¬ 
quired  to  obtain  this  central  deflection  may  be  sought. 

The  numerical  solution  may  be  more  lengthy  with  some  choices  of 
parameters  than  with  others.  If  p,  A0,  and  e0  are  given,  the  procedure  is 
simpler  than,  for  example,  when  p,  /i0 ,  and  a  are  given,  for  t0  can  be  comput¬ 
ed  from  e0  and  h0  by  Equation  [25]  and  each  integration  of  Equations  [26] 
gives  a  better  approximation  of  tr  and  e0 .  When  tr  and  ee  have  been  deter¬ 
mined  with  a  sufficient  degree  of  accuracy  as  functions  of  r,  the  radius  of 
the  plate  is  determined  by  setting  ee  equal  to  zero  and  solving  for  r.  This 
value  of  r  is  the  radius  a.'  The  deflection  z  is  then  found  by  integration 
using  Equation  [21],  and  z0  is  determined  so  that  z(a)  *  0.  If,  instead,  p, 

A0,  and  the  radius  a  of  the  plate  are  considered  known,  estimates  are  first 
made  for  z0,  e0,  t0.*  Then  Equations  [26],  [25],  and  [11]  are  used  itera¬ 
tively  to  find  sufficiently  accurate  determinations  of  tr  and  e0l  and  a  calcu¬ 
lated  radius  of  the  plate  a'  is  obtained  by  setting  ee  equal  to  zero.  On  the 
basis  of  these  results  new  evaluations  of  z0,  e0,  and  fc0  are  made  and  the 
process  is  repeated  until  a'  is  sufficiently  close  to  the  given  radius  a  of 
the  plate.  In  the  case  where  a  family  of  solution  curves  Is  desired,  it  is 
equally  effective  and  more  direct  to  proceed  by  the  first  method,  selecting 
suitable  sets  of  values  for  p,  h0,  and  e0. 

FIRST  APPROXIMATIONS 

To  facilitate  the  choice  of  first  approximations  for  eg  ,  tr,  etc., 
the  following  formulas  are  set  down.  Experiments  show  that  a  deflected  plate 
is  approximately  a  spherical  cap.  Consequently  the  tensions  t9  and  tT  are 
approximately  equal  and  constant  throughout  the  plate.  That  is,  the  relation¬ 
ship 

tr  ~  t$  —  to 

where  t0  is  the  common  tension  at  the  center,  is  approximately  true  for  all 
radial  values  r.  It  then  follows  from  Equation  [23]  that 

=  e6 

A  first-order  approximation  to  a  spherical  deflection  function  is 
the  parabolic  deflection  function  vanishing  at  r  =  a: 

z  -  J„(l  -  £)  127) 


#  Estimates  may  be  made  from  formulas  given  in  the  next  section. 
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where  x0  is  the  deflection  at  r  *  0 

It  may  be  shown  ( 1 4 )  that,  for  a  membrane  of  tension  t0>  p  is  given 
precisely  by 


V  = 


44*0 
a2  +  z0* 


This  form  is  preferable  when  z0  is  large  compared  to  a. 

Experiments  show  that  ee  varies  in  an  approximately  parabolic  man¬ 
ner.  Hence,  as  a  first  approximation 


where  c0  is  the  strain  at  the  center. 

The  quantity  t0/h0  can  be  calculated  when  e0  is  pre-assigned.  For 
at  r  *  0,  er  ~  e9  =  e0,  and  Equation  [25]  reduces  to 

-jk  =  6^  €0d  “  2e0) 

"'O  yo 


where  y0  is  the  value  of  y  at  r  =  0.  By  Equation  [11] 

y0  = 

and 


Yo  =  W  T(>'o)(1  “  2eo) 

Finally,  there  may  be  written  the  approximate  law 


Co  = 


The  plausibility  of  this  relationship  may  be  demonstrated  as  follows.  Com¬ 
bining  Equation  [ 1 6 ]  and  the  approximate  equations  er  =  *  e0(1  -  r2/a 2) 

which  were  found  in  the  foregoing,  we  obtain 


h  o  —  h  — 


i.e.,  a  parabolic  variation  of  h0  -  h  is  implied.  The  volume  V  of  the  plate 
before  deformation  is  equal  to  its  volume  after  deformation,  Hence, 


V 


*  ira2h0 
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Therefore,  to  a  first  approximation 

+  !(£)’][!- 2<o(l  -  £)]*-.• 


Since 


and 


it  follows  that 


a 

/[ 


2  -  *i 


i1  ~  h) 


(dzf  _  4r2z02 
\dri  a4 


2  r  + 


4  rV1 


1  -  2«„(l  -  £)' 


dr  =“  a2 


or 


/{2r  +  ^  -  4‘»r(1  -  $  -  ^[2'.(i  -  -pf*'-”2 

4^2  r  /  ^,2  n 

We  may  neglect  the  term  — ^T"r€o\1  ~  ^zjl  since  it  is  small  relative  to  the 
others.  The  equation  may  then  be  reduced  to 


f  f  4r3z02  .  /,  r2  \' 

J  rV^  -  M1  _  7)J 


dr  —  0 


Finally, 


«.*  -  mi  -  f) 


or 


=  K 

a2 


€0  -  ~2 


For  convenience  the  foregoing  approximate  formulas  are  rewritten  in 
collected  form 


*  - 

€0  _  n2 


a * 

4toZ0  4to^~ 

~ai  +  z02  *  ct(l  +  e0) 


[28] 
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h~0  "  TT  TW(1  ~  2€°} 


[28] 


[29] 


Equations  [28]  give  four  independent  relations  among  the  seven  quantities  t0, 
Zq  ,  e0,  h0,  p,  a,  and  y0.  Hence,  if  any  three  of  these  seven  quantities  are 
given,  the  others  may  be  determined  if  r(y)  is  known.  The  successive  approxi¬ 
mations  furnished  by  Equations  [26],  [25],  and  [11]  may  then  be  determined, 
starting  with  the  evaluations  for  tr  and  ee  based  on  Equations  [29]. 

This  completes  the  consideration  of  initial  approximations  for  use 
in  Equations  [26],  [25],  and  [11]. 


APPLICATION  OF  THE  METHOD 

The  methods  discussed  in  the  foregoing  for  predicting  the  action  of 
a  circular  clamped  plate  under  pressure  will  now  be  applied.  The  material 
selected  to  show  the  application  of  the  method  is  medium  steel  having  the 
r(y)  characteristic  of  Figure  3  for  the  tensile  specimen  with  the  stress- 
strain  curve  of  Figure  2;  see  pages  6  and  7  respectively.  For  the  purposes 
of  starting  the  computations  a  nominal  radius  of  10  inches  and  a  thickness  of 
1/8  inch  were  chosen  and  used  in  the  initial  formulas.  The  results  as  finally 
presented  in  Figures  8  to  12  inclusive  are  in  terms  of  the  general  parameters 
discussed  on  page  15,  and  so  apply  to  any  thin  plate  with  a  similar  v(y) 
characteristic . 

The  first  step  in  the  computation  is  to  select  a  range  of  values 
of  z 0  and  substitute  them  and  the  chosen  values  of  h0 ,  a,  and  r(y)  in  the 
approximate  formulas,  Equations  [28],  to  obtain  a  table  such  as  Table  1. 

The  second  step  in  the  computation  is  to  tabulate  tr  and  e6  for 
selected  values  of  radial  distance,  using  the  first  and  third  of  Equations 
[293.  The  values  of  e0*  *o»  anci  V  in  Table  1  are  now  to  be  considered  exact 
and  will  be  used  throughout  the  remaining  computations.  The  values  of  a  and 
*0  used  thus  far  are  entirely  discarded  from  this  point  on.  They  served 
merely  as  a  convenience  for  starting  the  calculations,  New  and  exact  values 
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TABLE  1 


Approximate  Values  of  «0,  t0,  and  p 

The  radius  a  is  taken  to  be  10  inches  and  h0  to  be  l/8  inch. 


*0 

in  inches 

,  _  £0! 
€°  a2 

<0  =  72  ho  T(y0)(l  “  2*o) 
in  pounds  per  inch 

JJ  _  ,4  Vzo 

p  a*  +  Zq* 

in  pounds  per 
square  inch 

0.2 

0.0004 

3206 

25.64 

0.5 

0.0025 

4670 

93.17 

1 .0 

0.01 

6289 

251 .56 

1.5 

0.0225 

7510 

4i6.5 

2.0 

0.04 

8758 

699.00 

2.5 

0.0625 

9354 

878.5 

5.0 

0.09 

9528 

1049.00 

of  z0  and  a  will  later  be  determined  for  each  set  of  values  e0,  t0,  and  p. 

The  third  step  is  the  graphical  determination  from  Figure  7  of  er  correspond¬ 
ing  to  the  values  of  tr  and  et  tabulated  in  Step  2.  The  fourth  step  is  to 
substitute  these  evaluations  of  tr,  ee,  and  er  in  the  integration  formulas 
[26]  to  secure  new  values  of  tr  and  e9  by  numerical  integration.  Steps  3  and 
4  are  then  repeated  until  successive  approximations  of  tr,  e9,  and  er  show 
little  change.  If  now  the  values  of  ce  for  a  set  of  values  e0,  t0,  hQ,  and 
p  from  Table  1  are  plotted  against  radial  distance,,  it  will  be  found  that  e9 
becomes  zero  at  a  radial  distance  which  is  found  to  approximate  10  inches, 
the  value  used  to  obtain  Table  1 .  The  radial  distance  at  which  e9  is  zero 
must  be  considered  the  exact  value  of  the  radius  a  corresponding  to  the  set 
<0,  t0  ,  and  p.  This  value  of  a  is  used,  of  course,  in  computing  quantities 
such  as  r/a  and  pa/h0  for  this  set.  Similarly  other  values  of  a  approximat¬ 
ing  10  inches  will  be  found  for  other  sets  of  e0,  t0  ,  and  p.  The  quantities 
€r  and  were  computed  in  this  manner  and  the  results  are  plotted  for  the 
various  sets  as  shown  in  Figures  8  and  9 . 

Other  quantities  such  as  t9,  art  and  a9  may  readily  be  calculated 
by  means  of  Equations  [23],  [16],  and  [19]*  The  quantities  or  and  ae  were 
determined  in  this  manner  and  the  results  are  shown  in  Figures  10  and  11 . 

The  deflection  z  at  radial  distance  r  may  be  computed  by  Equation 
[21],  The  values  of  z0  thus  obtained,  which  differ  somewhat  from  those  used 
in  Table  1,  are  the  exact  values  of  central  deflection  corresponding  to  the 


(Text  continued  on  page  27) 
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Figure  8  -  Radial  Strain  Plotted  against  Radial  Distance 

Radial  strains  were  computed  by  integration  of  equilibrium  conditions,  strain-displacement 
relations,  and  stress-strain  equations  of  plasticity.  The  load-strain  curve  of  Figure  2, 
taken  from  a  tensile  test  on  a  specimen  of  the  same  material  as  the  plate,  medium  steel, 
was  employed. 

Experimental  points  numbered  1  to  7  are  to  be  compared  with  points  on  the  theoretical  curves 
1  to  7  respectively. 

It  should  be  noted  that  in  this  set  of  curves  t0/a  was  used  as  the  parameter  instead  of 

pa/h„. 
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Figure  9  -  Circumferential  Strain  Plotted  against  Radial  Distance 

Circumferential  strains  were  computed  by  integration  of  equilibrium  conditions,  strain- 
displacement  relations,  and  stress-strain  equations  of  plasticity.  The  load-strain  curve 
of  Figure  2  obtained  from  a  tensile  test  on  a  specimen  taken  from  the  same  material  as  the 
plate,  medium  steel,  was  employed;  see  page  6. 

Experimental  points  numbered  1  to  7  are  to  be  compared  with  points  on  the  theoretical 
curve  1  to  7,  respectively. 

It  should  be  noted  that  in  this  set  of  curves  x0/a  was  used  as  the  parameter  instead  of 

p»Ao< 
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Figure  10  -  Radial  Stress  Plotted  against  Radial  Distance 

Radial  stresses  in  a  thin  circular  plate  were  computed  by  integration  of  equilibrium,  strain- 
displacement,  and  plasticity  laws.  The  load-strain  curve  of  Figure  2  obtained  from  a  tensile 
test  on  a  specimen  taken  from  the  same  material  as  the  experimental  plate,  medium  steel,  was 
employed. 

Experimental  points  numbered  3  to  7  are  to  be  compared  with  points  on  the  theoretical  curves 
3  to  7,  respectively. 


Circumferential  Stress  <r,  in  pounds  per  square  inch 
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Figure  11  -  Circumferential  Stress  Plotted  on  Radial  Distance 

Circumferential  stresses  were  computed  for  a  thin  circular  plate  by  integration  of 
equilibrium,  strain-displacement,  and  plasticity  laws.  The  load-strain  curve  of  Figure  2 
obtained  from  a  tensile  test  on  a  specimen  taken  from  the  same  material  as  the  experi¬ 
mental  plate,  medium  steel,  was  employed;  see  page  6. 

Experimental  points  numbered  3  to  7  are  to  be  compared  with  points  on  the  theoretical 
curves  3  to  7,  respectively. 
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too  high;  see  Figure  3  on  page  7.  biaxlal  Stresses  about  8  Per  cent 


seta  «0>  «0»  and  p.  Figure  12  shows  curves  of  deflection  plotted  against 
radial  distance  for  various  pressures  thus  obtained. 

Figures  8,  9,  10,  11,  and  12  are  plotted  by  use  of  the  proportion¬ 
ate  coordinates  introduced  and  established  by  the  equations  on  page  15.  The 
figures  therefore  apply  to  any  thin  circular  plate  of/ material  with  the  r(r) 
characteristic  of  Figure  3*  as  given  on  page  7. 

COMPARISON  OF  THEORY  AND  EXPERIMENT 

With  each  of  the  graphs  for  «r,  eg,  ari  ag,  and  z  computed  as  in 
the  foregoing  are  shown  corresponding  experimental  points.  The  data  from 
which  these  experimental  points  were  taken  were  obtained  at  the  Taylor  Model 


28 


Basin  (15)  on  a  circular  medium-steel  plate  nominally  20  inches  in  diameter 
and  1/8  inch  thick.  This  plate  was  machined  from  the  same  stock  as  the  ten¬ 
sile  specimen  which  yielded  the  r(r)  curve  used  in  the  calculations  of  this 
report.  The  fact  that  in  the  tests  the  rim  bounding  the  plate  had  a  curved 
edge  of  1/2-lnch  radius,  as  in  Figure  1,  made  necessary  a  correction  to  the 
plate  radius.  In  the  computations  of  the  experimental  data  a  value  of  10.25 
Inches  was  assumed  as  the  radius  a. 

As  the  test  plate  was  deformed  under  pressure,  measurements  were 
made  of  the  deflection  z  and  the  strains  ee  and  er  at  selected  radii.  From 
these  data  and  Equations  [16],  [17],  [  1 8 ] ,  and  [19]  the  values  of  <rr  and  a9 
were  computed  and  plotted  as  the  experimentally  determined  points. 

The  deflection  curves  of  Figures  12  and  15  are  considered  to  be  the 
most  significant  check  of  theory  and  experiment  in  that  the  deflection  is  sub¬ 
ject  to  direct  observation  with  only  small  error.  The  two  sets  of  curves  in 
Figures  12  and  15  correspond  in  pairs  to  selected  values  of  the  proportionate 
pressure  pa/h0 .  It  will  be  seen  that  at  low  pressures  the  agreement  is  within 
a  few  per  cent.  As  the  pressure  is  increased,  the  experimental  values  exceed 


Figure  15  -  Comparison  of  Theoretical  and  Experimental  Curves 
of  Pressure  against  Center  Deflection 

The  theoretical  curve  was  obtained  by  integration  of  equilibrium  strain-displacement, 
and  plasticity  laws  using  a  load-strain  curve  obtained  from  a  medium-steel  specimen. 

The  experimental  curve  was  obtained  by  direct  measurement  on  a  medium-steel  plate  of 
nominal  radius  10  inches  and  thickness  1/8  inch. 
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the  calculated  values  by  larger  amounts;  the  largest  discrepancy  is  about  8 
per  cent. 

The  lower  of  the  two  r(y)  curves  in  Figure  3  on  page  7  was  obtained 
from  the  experimental  stresses  and  strains  for  the  center  of  the  plate,  by 
the  use  of  Equations  [2]  and  [11].  This  curve  may  be  referred  to  as  the 
experimental  r{y)  curve  for  the  plate.  The  other  curve  shown  in  Figure  3, 
although  obtained  experimentally  from  a  uniaxial  test,  may  conveniently  be 
regarded  as  the  theoretical  r(y)  curve  for  the  plate.  It  will  be  noted  that 
the  theoretical  values  of  r  exceed  the  experimental  ones  for  a  given  value 
of  y.  Because  of  high  theoretical  values  of  r,  theoretical  stresses  should 
be  high  for  a  given  pair  of  strains.  This  statement  is  corroborated  by  the 
stress  values  shown  in  Figures  10  and  11.  As  expected,  the  theoretical  val¬ 
ues  of  stress  are  generally  higher  than  the  experimental  ones  for  a  given 
value  of  pa/h0;  at  high  pressures  the  difference  may  be  as  much  as  10  per 
cent.  As  a  result  of  high  theoretical  stresses,  the  theoretical  pressure 
is  high  for  a  given  deflection.  It  follows  that  the  theoretical  deflections 
should  be  lower  than  the  experimental  ones  for  a  given  pressure.  These  de¬ 
ductions  are  also  corroborated  by  the  deflection  curves  of  Figures  12  and  13. 

In  Figures  8  and  9  on  pages  23  and  24  the  theoretical  and  experi¬ 
mental  values  of  strain  are  presented  for  selected  values  of  z0/a.  It  is 
believed  that  a  major  portion  of  the  discrepancy  between  theoretical  and  ex¬ 
perimental  values  of  strain  may  be  ascribed  to  errors  arising  from  difficul¬ 
ties  encountered  in  using  metalectric  strain  gages.  It  will  be  seen  that  at 
low  deflections  the  theoretical  values  of  strain  usually  are  exceeded  by  the 
computed  values.  The  measured  values  exceed  estimates  based  on  any  of  the 
known  plate  or  membrane  formulas  for  small  or  moderate  deflections.  As  the 
deflection  is  increased,  there  is  agreement  on  tne  average  in  Figure  8  ex¬ 
cept  at  large  values  of  radial  distance.  A  small  portion  of  the  discrepancy 
near  the  edge  of  the  plate  may  be  due  to  bending  strains  since  such  strains 
were  not  taken  into  account  ip  the  theory.  In  this  region  the  experimental 
plate  should  be  less  convex  than  the  theory  indicates,  and  the  measured  ra¬ 
dial  strains  accordingly  less  for  a  given  value  of  z0/a  . 

Computations  not  described  in  this  report  show  that  a  better  check 
of  theoretical  and  experimental  values  of  deflection  will  result  if  logarith¬ 
mic  strains  are  used  in  place  of  the  conventional  strains  in  the  assumed  plas¬ 
ticity  laws.  Such  a  refinement  and  the  elimination  of  the  discrepancy  between 
uniaxial  and  biaxial  r(y)  characteristics  obtained  in  the  tests  are  believed 
generally  to  outweigh  other  refinements  that  could  be  made  in  the  theory, 
such  as  the  inclusion  of  bending  stresses  or  the  use  of  exact  strain- 
displacement  relations  in  the  computations  in  place  of  first-order  approx¬ 
imations  . 
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CONCLUSIONS 

The  plasticity  laws,  Equations  [1]  and  [2],  together  with  equilib¬ 
rium  conditions  and  strain-displacement  relationships  for  the  clamped,  thin, 
circular  plate  under  pressure,  may  be  numerically  integrated  to  give  theoreti¬ 
cal  values  of  stresses,  strains,  and  deflections  at  a  given  pressure. 

Results  so  obtained  were  found  to  be  in  agreement  with  experiment 
within  about  10  per  cent. 
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